**Civil Military Relations, Postwar Politics, and Civil War Recurrence
**Louis-Alexandre Berg, Georgia State University, LBerg@gsu.edu


use "PSGD_Civil War Recurrence Replication Data_Final" 

*Figure 1
catplot officer_appointment if finalyear==1, percent stack asyvars
graph save Graph ApptGraph.gph 

catplot civ_access if finalyear==1, percent stack asyvars
graph save Graph AccessGraph.gph 

catplot civ_control  if finalyear==1, percent stack asyvars
graph save Graph civ_controlGraph.gph 

gr combine ApptGraph.gph  AccessGraph.gph   civ_controlGraph.gph   

***Table A.3: Summary Statistics
sum officer_appointment   civ_access civ_control mil_integration milint_full milexppercapita milexppcchng_oneyear lngdpcap  lnpop  ethfrac resourcerents victory  peaceagmt duration  factions un polity2 coldwar 

***Cox Hazard Model, by peacespell

**Table 1: Main Analysis
stset peaceyear, failure(recur_high) id(peacespell)
stcox   officer_appointment , cluster(id) robust
estimates store Model_1
stcox   officer_appointment lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar , cluster(id) robust
estimates store Model_2
stcox   civ_access, cluster(id) robust
estimates store Model_3
stcox   civ_access lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  , cluster(id) robust
estimates store Model_4
stcox   civ_control  , cluster(id) robust
estimates store Model_5
stcox   civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  , cluster(id) robust
estimates store Model_6
stcox   civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions    coldwar  , cluster(id) robust
estimates store Model_7
stcox   officer_appointment civ_access civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions  polity2  coldwar  , cluster(id) robust
estimates store Model_8

esttab     ///
   Model_1  ///
   Model_2  ///
   Model_3  ///
   Model_4  ///
   Model_5  ///
  Model_6 ///
   Model_7 ///
   Model_8 ///
   using table2.rtf, ///
 eform se(3)  starlevels(* 0.10 ** 0.05 *** 0.01) stats(N  N_sub N_clust N_fail ll  ) compress nogap  replace label   nolabel order(officer_appointment civ_access civ_control  lngdpcap  lnpop  ethfrac  resourcerents    peaceagmt victory  un duration factions  polity2 coldwar)
 
 **Figure 2
 stset peaceyear, failure(recur_high) id(peacespell)
stcox   officer_appointment lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar , cluster(id) robust
stcurve, survival at1(officer_appointment=1) at2(officer_appointment=0) graphregion(lcolor(white) ilcolor(white)   ifcolor(white) )
graph save  figure2_1.gph 
stcox   civ_access lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  , cluster(id) robust
stcurve, survival at1(civ_access=1) at2(civ_access=0)
graph save  figure2_2.gph 
stcox   civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  , cluster(id) robust
stcurve, survival at1(civ_control =1) at2(civ_control =0)
graph save  figure2_3.gph 

**Survival Curves with Confidence Intervals, no controls, included in Supplementary Appendix Figure A.1
sts graph, survival ci by(officer_appointment) name(g1, replace) title ("Officer officer_appointment") graphregion(lcolor(white) ilcolor(white)   ifcolor(white) )
 sts graph, survival ci by(civ_access) name(g2, replace) title ("Civilian Access") graphregion(lcolor(white) ilcolor(white)   ifcolor(white) )
sts graph, survival ci by(civ_control ) name(g3, replace) title ("Civilian Control") graphregion(lcolor(white) ilcolor(white)   ifcolor(white) )
  

**TABLE 3: INTEGRATION AND CAPACITY
 estimates drop Model_9 Model_10 Model_11 Model_12 Model_13 Model_14   
 stset peaceyear, failure(recur_high) id(peacespell)
stcox  mil_integration lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar , cluster(id) robust
estimates store Model_9
stcox   milint_full lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  , cluster(id) robust
estimates store Model_10
stcox    milint_full officer_appointment civ_access civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar , cluster(id) robust
estimates store Model_11
stcox   milexppercapita lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  , cluster(id) robust
estimates store Model_12
stcox   milexppcchng_oneyear lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  , cluster(id) robust
estimates store Model_13
stcox   milexppercapita milexppcchng_oneyear officer_appointment FormalOversight civ_control   lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar , cluster(id) robust
estimates store Model_14

esttab     ///
   Model_9  ///
   Model_10  ///
   Model_11 ///
   Model_12  ///
   Model_13  ///
  Model_14 ///
  using table3.rtf, ///
 eform se(3)  starlevels(* 0.10 ** 0.05 *** 0.01) stats(N  N_sub N_clust N_fail ll  ) compress nogap  replace label   nolabel order(mil_integration milint_full milexppercapita   milexppcchng_oneyear  officer_appointment civ_access civ_control  FormalOversight   lngdpcap  lnpop  ethfrac  resourcerents    peaceagmt victory  un duration factions  polity2 coldwar)
  
 **Figure 3

 gen officerfile = .
replace  officerfile = 1 if officer_appointment==1 & milint_full==1
replace  officerfile = 0 if officer_appointment==0
replace  officerfile = 0 if milint_full==0

sts graph, survival ci by(officer_appointment) name(g4, replace) title ("Officer officer_appointment") graphregion(lcolor(white) ilcolor(white)   ifcolor(white) )
graph save  figure3_1.gph 
sts graph, survival ci by(milint_full) name(g5, replace) title ("Rank and File Integration")  graphregion(lcolor(white) ilcolor(white)   ifcolor(white) )
graph save  figure3_2.gph
sts graph, survival ci by(officerfile) name(g6, replace) title ("Officer officer_appointment & Rank and File Integration")  graphregion(lcolor(white) ilcolor(white)   ifcolor(white) )
graph save  figure3_3.gph
 
 ***ROBUSTNESS CHECKS
 
***Table A.4: Weibull Model 
estimates drop Model_1 Model_2 Model_3 Model_4 Model_5 Model_6 Model_7 Model_8
stset peaceyear, failure(recur_high) id(peacespell)
streg   officer_appointment , distribution (weibull) cluster(id) robust
estimates store Model_1
streg   officer_appointment lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar , distribution (weibull) cluster (id) robust
estimates store Model_2
streg   civ_access if  oneyear !=1,distribution (weibull)   cluster(id) robust
estimates store Model_3
streg   civ_access lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  ,distribution (weibull)  cluster(id) robust
estimates store Model_4
streg   civ_control  , distribution (weibull) cluster (id) robust
estimates store Model_5
streg   civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  , distribution (weibull) cluster (id) robust
estimates store Model_6
streg   civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions    coldwar  , distribution (weibull) cluster(id) robust
estimates store Model_7
streg   officer_appointment civ_access civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions    coldwar  , distribution (weibull) cluster (id) robust
estimates store Model_8
esttab      ///
   Model_1  ///
   Model_2  ///
   Model_3  ///
   Model_4  ///
   Model_5  ///
  Model_6 ///
   Model_7 ///
   Model_8 ///
   using table4.rtf, ///
 eform se(3)  starlevels(* 0.10 ** 0.05 *** 0.01) scalars(N_sub N_clust) compress nogap  replace label   nolabel order(officer_appointment civ_access civ_control  lngdpcap  lnpop  ethfrac  resourcerents    peaceagmt victory un duration   factions   polity2 coldwar)

 
**Table A.5: Logit Model: Conflict Recurrence
estimates drop Model_1 Model_2 Model_3 Model_4 Model_5 Model_6 Model_7 Model_8

logit recur_high    officer_appointment time time2 time3 ,  cluster(id) robust
estimates store Model_1
logit recur_high     officer_appointment lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  time time2 time3 ,  cluster (id) robust
estimates store Model_2
logit recur_high     civ_access time time2 time3 if  oneyear !=1,  cluster(id) robust
estimates store Model_3
logit recur_high     civ_access lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar time time2 time3  , cluster(id) robust
estimates store Model_4
logit recur_high     civ_control  time time2 time3 ,  cluster (id) robust
estimates store Model_5
logit recur_high     civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  time time2 time3 ,  cluster (id) robust
estimates store Model_6
logit recur_high      civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions    coldwar time time2 time3  ,cluster(id) robust
estimates store Model_7
logit recur_high      officer_appointment civ_access civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions    coldwar time time2 time3  ,  cluster (id) robust
estimates store Model_8

esttab     ///
   Model_1  ///
   Model_2  ///
   Model_3  ///
   Model_4  ///
   Model_5  ///
  Model_6 ///
   Model_7 ///
   Model_8 ///
   using tableA5.rtf, ///
  se(3)  starlevels(* 0.10 ** 0.05 *** 0.01) stats(N  N_sub N_clust N_fail ll  ) compress nogap  replace label   nolabel order(officer_appointment civ_access civ_control  lngdpcap  lnpop  ethfrac  resourcerents    peaceagmt victory  un duration factions  polity2 coldwar)

**Table A.6: Cox Proportional Hazard Model with alternative specification of independent variables
 estimates drop Model_1 Model_2 Model_3 Model_4 Model_5 Model_6 Model_7 Model_8
stset peaceyear, failure(recur_high) id(peacespell)
stcox   officer_appointment , cluster(id) robust
estimates store Model_1
stcox   officer_appointment lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar , cluster(id) robust
estimates store Model_2
stcox   FormalOversight if  oneyear !=1, cluster(id) robust
estimates store Model_3
stcox   FormalOversight lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  , cluster(id) robust
estimates store Model_4
stcox   EffectiveOversight , cluster(id) robust
estimates store Model_5
stcox   EffectiveOversight lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  , cluster(id) robust
estimates store Model_6
stcox   EffectiveOversight lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions    coldwar  , cluster(id) robust
estimates store Model_7
stcox   officer_appointment FormalOversight EffectiveOversight lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions  polity2  coldwar  , cluster(id) robust
estimates store Model_8

esttab     ///
   Model_1  ///
   Model_2  ///
   Model_3  ///
   Model_4  ///
   Model_5  ///
  Model_6 ///
   Model_7 ///
   Model_8 ///
   using tableA6.rtf, ///
 eform se(3)  starlevels(* 0.10 ** 0.05 *** 0.01) stats(N  N_sub N_clust N_fail ll  ) compress nogap  replace label   nolabel order(officer_appointment FormalOversight EffectiveOversight  lngdpcap  lnpop  ethfrac  resourcerents    peaceagmt victory  un duration factions  polity2 coldwar)

 
 
 **Table A.7: Hazard Model with Additional Controls: Conflict Recurrence
 
 estimates drop Model_1 Model_2 Model_3 Model_4

stset peaceyear, failure(recur_high) id(peacespell)
stcox   officer_appointment exclpop  oilexp lngdpcap  lnpop      duration peaceagmt victory   un factions   polity2 coldwar , cluster(id) robust
estimates store Model_1
stcox   civ_access exclpop  oilexp lngdpcap  lnpop   duration peaceagmt victory   un factions   polity2 coldwar  , cluster(id) robust
estimates store Model_2
stcox   civ_control  exclpop  oilexp lngdpcap  lnpop       duration peaceagmt victory   un factions   polity2 coldwar  , cluster(id) robust
estimates store Model_3
stcox   officer_appointment civ_access civ_control  exclpop  oilexp lngdpcap  lnpop       duration peaceagmt victory   un factions   polity2 coldwar  , cluster(id) robust
estimates store Model_4

esttab     ///
   Model_1  ///
   Model_2  ///
   Model_3  ///
   Model_4 ///
   using tableA7.rtf, ///
 eform se(3)  starlevels(* 0.10 ** 0.05 *** 0.01) stats(N  N_sub N_clust N_fail ll  ) compress nogap  replace label   nolabel order(officer_appointment civ_access civ_control  lngdpcap  lnpop  exclpop  oilexp   peaceagmt victory  un duration factions  polity2 coldwar)

 
***Table A.8: Random Effects Model: Conflict Recurrence
  estimates drop Model_1 Model_2 Model_3 Model_4 Model_5 Model_6 Model_7 Model_8
xtset id year
xtlogit recur_high    officer_appointment  ,  re
estimates store Model_1
xtlogit recur_high     officer_appointment lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar , re
estimates store Model_2
xtlogit recur_high     civ_access  if  oneyear !=1,  re
estimates store Model_3
xtlogit recur_high     civ_access lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  , re
estimates store Model_4
xtlogit recur_high     civ_control   ,  re
estimates store Model_5
xtlogit recur_high     civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  ,  re
estimates store Model_6
xtlogit recur_high      civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions    coldwar  , re
estimates store Model_7
xtlogit recur_high      officer_appointment civ_access civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions    coldwar  ,  re
estimates store Model_8
esttab      ///
   Model_1  ///
   Model_2  ///
   Model_3  ///
   Model_4  ///
   Model_5  ///
  Model_6 ///
   Model_7 ///
   Model_8 ///
   using tableA8.rtf, ///
    se(3)  starlevels(* 0.10 ** 0.05 *** 0.01) scalars(N_sub N_clust) compress nogap  replace label   nolabel order(officer_appointment civ_access civ_control  lngdpcap  lnpop  ethfrac  resourcerents  peaceagmt victory  un duration  factions   polity2 coldwar)
 
***Table A.9: Cox Proportional Hazard Model with Region Fixed Effects
 
 estimates drop Model_1 Model_2 Model_3 Model_4 Model_5 Model_6 Model_7 Model_8

stset peaceyear, failure(recur_high) id(peacespell)
stcox   officer_appointment i.Region  , cluster(id) robust
estimates store Model_1
stcox   officer_appointment lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar i.Region  , cluster(id) robust
estimates store Model_2
stcox   civ_access i.Region  if  oneyear !=1, cluster(id) robust
estimates store Model_3
stcox   civ_access lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar i.Region   , cluster(id) robust
estimates store Model_4
stcox   civ_control  i.Region , cluster(id) robust
estimates store Model_5
stcox   civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar i.Region    , cluster(id) robust
estimates store Model_6
stcox   civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions    coldwar i.Region   , cluster(id) robust
estimates store Model_7
stcox   officer_appointment civ_access civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions  polity2  coldwar i.Region   , cluster(id) robust
estimates store Model_8

esttab      ///
   Model_1  ///
   Model_2  ///
   Model_3  ///
   Model_4  ///
   Model_5  ///
  Model_6 ///
   Model_7 ///
   Model_8 ///
   using tableA9.rtf, ///
 eform se(3)  starlevels(* 0.10 ** 0.05 *** 0.01) scalars(N_sub N_clust) compress nogap  replace label   nolabel order(officer_appointment civ_access civ_control  lngdpcap  lnpop  ethfrac  resourcerents    peaceagmt victory un duration  factions   polity2 coldwar)
 


***Table A.10: Cox Proportional Hazard Model with Decade Fixed Effects: Conflict Recurrence

estimates drop Model_1 Model_2 Model_3 Model_4 Model_5 Model_6 Model_7 Model_8
stset peaceyear, failure(recur_high) id(peacespell)
stcox   officer_appointment  sixties seventies eighties nineties oughts , cluster(id) robust
estimates store Model_1
stcox   officer_appointment lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar   sixties seventies eighties nineties oughts , cluster(id) robust
estimates store Model_2
stcox   civ_access   sixties seventies eighties nineties oughts, cluster(id) robust
estimates store Model_3
stcox   civ_access lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  sixties seventies eighties nineties oughts   , cluster(id) robust
estimates store Model_4
stcox   civ_control    sixties seventies eighties nineties oughts , cluster(id) robust
estimates store Model_5
stcox   civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  sixties seventies eighties nineties oughts  , cluster(id) robust
estimates store Model_6
stcox   civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions    coldwar   sixties seventies eighties nineties oughts  , cluster(id) robust
estimates store Model_7
stcox   officer_appointment civ_access civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions  polity2  coldwar   sixties seventies eighties nineties oughts  , cluster(id) robust
estimates store Model_8

esttab      ///
   Model_1  ///
   Model_2  ///
   Model_3  ///
   Model_4  ///
   Model_5  ///
  Model_6 ///
   Model_7 ///
   Model_8 ///
   using tableA10.rtf, ///
  eform se(3)  starlevels(* 0.10 ** 0.05 *** 0.01) scalars(N_sub N_clust) compress nogap  replace label   nolabel order(officer_appointment civ_access civ_control  lngdpcap  lnpop  ethfrac  resourcerents    peaceagmt victory un duration  factions   polity2 coldwar sixties seventies eighties nineties oughts )
 
 
***Table A.11 & Figure A.3: Logit Model with interactions with time: Civil War Recurrence
 estimates drop Model_1 Model_2 Model_3 Model_4 Model_5 Model_6 Model_7 Model_8

logit recur_high    officer_appointment##c.peaceyear time time2 time3 ,  cluster(id) robust
estimates store Model_1
margins, dydx(officer_appointment) at(peaceyear = (0(1)50))
marginsplot
logit recur_high     civ_access##c.peaceyear time time2 time3,  cluster(id) robust
estimates store Model_2
margins, dydx(civ_access) at(peaceyear = (0(1)50))
marginsplot
logit recur_high     civ_control##c.peaceyear time time2 time3 ,  cluster (id) robust
estimates store Model_3
margins, dydx(civ_control ) at(peaceyear = (0(1)50))
marginsplot


esttab     ///
   Model_1  ///
   Model_2  ///
   Model_3  ///
   using tableA11.rtf, ///
   se(3)  starlevels(* 0.10 ** 0.05 *** 0.01) stats(N  N_sub N_clust N_fail ll  ) compress nogap  replace label   nolabel order(officer_appointment civ_access civ_control  lngdpcap  lnpop  ethfrac  resourcerents    peaceagmt victory  un duration factions  polity2 coldwar)



***Table A.12 Logit Model with Short vs. Long Peace Spells: Conflict Recurrence  
  
estimates drop Model_1 Model_2 Model_3

stset peaceyear, failure(recur_high) id(peacespell)
logit recur_high    officer_appointment time time2 time3 if shortspell==1 ,  cluster(id) robust
estimates store Model_1
logit recur_high    officer_appointment time time2 time3 if longspell==1 ,  cluster(id) robust
estimates store Model_2
logit recur_high    civ_access time time2 time3 if shortspell==1,  cluster(id) robust
estimates store Model_3
logit recur_high    civ_access time time2 time3 if longspell==1 ,  cluster(id) robust
estimates store Model_4
logit recur_high    civ_control  time time2 time3 if shortspell==1,  cluster(id) robust
estimates store Model_5
logit recur_high    civ_control  time time2 time3 if longspell==1,  cluster(id) robust
estimates store Model_6

esttab     ///
   Model_1  ///
   Model_2  ///
   Model_3  ///
   Model_4  ///
   Model_5  ///
  Model_6 ///
using tableA12.rtf, ///
 se(3) starlevels(* 0.10 ** 0.05 *** 0.01) stats(N  N_sub N_clust N_fail ll  ) compress nogap  replace label   nolabel order(officer_appointment civ_access civ_control  time time2 time3)

 
 
***Table A.13: Cox Proportional Hazard Model: Conflict Recurrence based on 25 battle death threshold
 
 estimates drop Model_1 Model_2 Model_3 Model_4 Model_5 Model_6 Model_7 Model_8

stset peaceyear_low, failure(recur_low) id(peacespell_low)
stcox   officer_appointment, cluster(id) robust
estimates store Model_1
stcox   officer_appointment lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar , cluster(id) robust
estimates store Model_2
stcox   civ_access , cluster(id) robust
estimates store Model_3
stcox   civ_access lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar  , cluster(id) robust
estimates store Model_4
stcox   civ_control  , cluster(id) robust
estimates store Model_5
stcox   civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions   polity2 coldwar , cluster(id) robust
estimates store Model_6
stcox   civ_control  lngdpcap  lnpop  ethfrac  resourcerents   duration peaceagmt victory   un factions    coldwar  , cluster(id) robust
estimates store Model_7
stcox   officer_appointment civ_access civ_control  lngdpcap  lnpop  ethfrac  resourcerents    peaceagmt victory   un duration factions  polity2  coldwar  , cluster(id) robust
estimates store Model_8
esttab      ///
   Model_1  ///
   Model_2  ///
   Model_3  ///
   Model_4  ///
   Model_5  ///
  Model_6 ///
   Model_7 ///
   Model_8 ///
  using table20.rtf, ///
 eform se(3)  starlevels(* 0.10 ** 0.05 *** 0.01) scalars(N_sub N_clust) compress nogap  replace label   nolabel order(officer_appointment civ_access civ_control  lngdpcap  lnpop  ethfrac  resourcerents    peaceagmt victory   un duration factions  polity2  coldwar)

 
 
